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Conversion  Factors,  Non-SI  to  SI  (Metric) 
Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  nan  be  converted  to  SI 
(metric)  units  as  follows; 


Multiply 


B Z 


To  Obtain 


feet 

inches 

inches  per  second  squared 

inch-kips  (force) 

inch-pounds  (force) 

kips  (1000  lb  force) 

pound  (force)-feet 
squared 

pound  (force) -inches 
squared 

pounds  (force) 

pound  (force) -seconds 
squared  per  inch 

pound  (force)-seconds 
squared  per  square  foot 

standard  free  fail  (g's) 


0.3048 

0.0254 

0.0254 

112.9848 

0.1129848 

4448.222 

0.413253 

0.028693 

4.448222 

175.1268 

47.88026 

9.80665 


mecres 

metres 

metres  per  second  squared 

newton-metres 

newton-metres 

newtons 

newton-metres  squared 
newton-metres  squared 
newtons 

newton-seconds  squared 
per  metre 

pascal-seconds  squared 
metres  per  second  squared 
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SEISMIC  ANALYSIS  OF  INTAKE  TOWERS 


Introduction 


1.  Intake  towers  in  the  outlet  works  of  many  existing  and 
planned  Corps  of  Engineers  projects  are  subjected  to  earthquake  hazards 
ranging  from  minor  to  severe.  The  Corps  has  traditionally  evaluated 
the  seismic  safety  of  these  intake  towers  through  a  dated,  gross  ap¬ 
proximation  of  earthquake  loading  known  as  the  seismic  coefficient 
method.  The  research  discussed  in  this  report  was  undertaken  to  de¬ 
velop  practical  procedures  fer  analysis  of  intake  tower  seismic  re¬ 
sponse  which  would  be  consistent  with  the  current  understanding  of 
earthquake  loading  and  structural  response.  However,  the  analytical 
procedures  developed  require  that  some  undesirable  approximations  of 
seismic  response  be  employed.  These  are  recommended  as  topics  for 
further  study. 


Seismic  Coefficient  Method 


2.  The  seismic  coefficient  method  (Office,  Chief  of  Engineers 
1964)  is  a  procedure  used  to  calculate  an  approximate  static  repre¬ 
sentation  of  earthquake  forces  which  can  be  applied  to  intake  towers 
and  other  structures  for  use  in  structural  analyses.  The  earthquake 
load  is  computed  by  multiplying  a  dimensionless  seismic  coefficient 
by  the  distributed  weight  of  the  intake  tower.  The  magnitude  of  the 
seismic  coefficient  is  determined  by  the  seismic  hazard  associated 
with  generalized  geographic  regions.  Once  the  static  representation  of 
the  earthquake  forces  has  been  calculated,  it  is  applied  to  the  intake 
tower,  and  the  shears  and  moments  are  computed  from  equilibrium 
requirements . 

3.  A  major  deficiency  of  the  seismic  coefficient  method  i:-  the 
inherent  assumption  of  rigid  body  motion.  Since  an  intake  tower  is 
not  rigid,  this  procedure  ignores  the  dependency  cf  structure  response 
on  the  natural  frequencies  of  the  tower  and  the  frequency  content  of 
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the  earthquake  gvoucd  motion.  Another  shortcoming  is  the  lack  of 
accuracy  in  tne  specification  of  the  seismic  coefficient.  The  same 
value  is  assigned  to  a  large  geographic  region.  Hence,  it  does  not 
account  for  the  seismic  characteristics  of  particular  locations. 
Therefore,  considering  these  limitations,  the  seismic  coefficient 
method  of  analysis  cannot  be  expected  to  accurately  predict  the  re¬ 
sponse  cf  an  intake  tower  subjected  to  an  earthquake  loading. 

Design  Earthquake 

4.  The  first  step  in  seismic  design  is  choosing  the  level  of  the 
seismic  load  effect  an^,  structure  might  experience  during  its  life. 

That  this  effect  cannot  be  adequately  represented  by  a  common  static 
load  for  all  structures  can  be  shown  by  examining  the  sing Le-degree- 
of-freedom  (SDOF)  idealization  of  a  structure  in  Figure  1.  If 
m  =  mass  of  the  structure 
c  =  damping  constant  of  the  structure 
k  =  stiffness  of  the  structure 

u(t)  =  displacement  of  the  structure  relative  to  the 
ground 

u(t)  =  velocity  of  the  structure  relative  to  the  ground 

ii(t)  =  acceleration  of  the  structure  relative  to  the 
ground 

ii  (t)  =  ground  acceleration 
the  equation  of  motion  of  this  structure  is 


mii  +  cu  +  ku  =  -mu  (t) 

8 


(1) 


Therefore,  the  internal  seismic  force  felt  by  the  structure 


ku  =  -m[ii  (t)  +  ii ( t )  ]  -  cu(t) 
8 


(2) 


is  a  time-dependent  function  of  its  stiffness,  damping,  and  inertial 
properties  as  well  as  of  the  earthquake  ground  motion. 


Figure  i.  Single-degree-of-f reedom 
idealization  of  a  structure 


j — u{t) 


5.  However,  Equation  1  does  suggest  that  one  method  of  repre¬ 
senting  the  seismic  force  is  through  the  specification  of  time  histo¬ 
ries  of  ground  acceleration  ii  (t)  .  The  time  histories  may  be  taken 

8 

from  actual  earthquakes,  such  as  shown  in  Figure  2,  or  may  be  randomly 
generated  by  computer  algorithms.  This  representation  of  the  earth¬ 
quake  loading  is  usually  necessary  for  structures  characterized  by  a 
nonlinear  response. 


Figure  2.  Accelerogram  from  the  El  Centro  earthquake, 

18  May  1940  (NS  component) 

6.  The  response  spectrum  is  another  way  of  specifying  the  seis¬ 
mic  load  effect.  A  response  spectrum  is  a  graphical  representation  of 
the  maximum  response  of  all  linear  SDOF  systems  for  a  specific  ground 
motion  at  a  particular  level  of  damping  (Clough  and  Penzien  1975). 
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For  each  time  history,  there  is  a  unique  response  spectrum  as  shown 
in  Figure  3.  A  response  spectrum  for  use  in  design  is  sometimes 
obtained  by  smoothing  the  spectra  of  several  actual  time  histories. 


Figure  3.  Response  spectra  for  the  El  Centro  earthquake,  18  May  1940 

(NS  component) 


7.  There  are  a  number  of  accepted  methods  at  present  for  deter¬ 
mining  design  earthquakes  for  different  location  (Hays  1980;  Kri- 
nitzsky*).  These  methods  consider  the  following  factors  in  prescribing 
an  ensemble  of  time  histories  and/or  the  response  spectrum  appropriate 
for  a  particular  site: 

a.  Geologic,  geophysical,  and  seismological  data  about  the 
site. 

b.  Presence  or  absence  of  recognizable  active  faults. 

c.  Estimated  upper  bound  magnitudes  for  earthquakes 
that  might  be  produced  by  these  faults. 

d.  Peak  motions  associated  with  these  events. 

e.  Selection  of  a  proper  attenuation  of  these  motions  from 
source  to  site. 

f.  Effects  of  site  characteristics  on  these  motions. 

g.  Likelihood  of  exceeding  various  ground  motions  using 
probabilistic  procedures. 

While  a  structural  engineer  is  primarily  interested  in  the  prescribed 
design  earthquake,  he  should  appreciate  the  professional  judgments 


*  E.  L.  Krinitzsky,  "Essentials  for  Specifying  Earthquake  Motions  in 
Engineering  Design,"  U.  S.  Army  Engineer  Waterways  Experiment  Sta¬ 
tion,  Vicksburg,  Miss.,  1981. 


involved  in  the  procedure  that  led  to  this  prescription.  Such  an 
appreciation  will  permit  him  to  appropriately  tailor  the  conservatism 
of  his  structural  analysis  to  assure  the  public  a  safe  and  economical 
project . 


Axisymmetric  Analyses 

8.  In  this  section,  analytical  models  of  axisymmetric  intake 
tower  response  to  earthquake  excitation  will  be  described  and  illus¬ 
trated.  It  will  be  seen  that  the  more  refined  analyses  may  require 
greater  effort  than  is  justified  in  many  cases.  Accordingly,  simpler 
but  less  accurate  methods  will  also  be  discussed  for  use  in  such 
situations . 

Analysis  with  program  EATSW 

9.  The  most  sophisticated  model  presently  available  considers 
an  axisymmetric  elastic  tower  surrounded  by  an  infinitely  extending 
reservoir  and  standing  on  a  rigid  foundation  (Liaw  and  Chopra  1975). 

A  special-purpose  finite  element  computer  code  (EATSW)  employing 
substructuring  and  fast  Fourier  transform  techniques  has  been  devel¬ 
oped  to  implement  this  analytical  model  (Liaw  and  Chopra  1973).  E/TSW 
numerically  evaluates  the  response  of  tins  analytical  model  to  the 
horizontal  component  of  an  earthquake  ground  motion  using  either 
explicit  mathematical  solutions  if  the  structure-water  interface  is 
cylindrical  or  a  finite  element  system  if  the  interface  geometry  is 
more  complex.  This  code  is  capable  of  including  any  number  modes 
of  vibration  desired  by  the  user.  EATSW  outputs  the  complete  time 
history  of  the  displacements  and  stresses  throughout  the  structure. 
Figures  4-13  show  two  different  intake  towers,  two  different  seismic 
motions,  and  the  responses  of  these  towers  to  these  motions  as  calcu¬ 
lated  in  Li.3w  and  Chopra  (1973).  In  the  EATSW  analyses,  only  three 
modes  of  vibration  of  the  San  Bernardino  Tower  and  four  modes  of  vib¬ 
ration  of  the  Briones  Tower  were  found  to  be  adequate.  For  subsequent 
comparison  with  simpler  analyse  ,  the  stress  output  from  EATSW  was 
converted  to  shears  and  moments  using  the  relationships 
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a.  SAN  BERNARDINO  TOWER  b.  FINITE  ELEMENT  IDEALIZATION 

Figure  4.  San  Bernardino  towtr  and  its  finite  element  idealization 
(elevations  (EL)  are  in  feet  referred  to  mean  sea  level) 


BRIONES  TOWER  b.  FINITE  ELEMENT  c.  FINITE  ELEMENT' 

IDEALIZATION  IDEALIZATIONS  THE 

OF  THE  TOWER  SURROUNDING  WATER 


Figure  5.  Briones  tower  and  finite  element  idealizations 
for  the  tower  and  surrounding  water 
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PSEUDOABSOLUTE  ACCELERATION  Sa,  g's 


TIME,  SEC 


*.  ACCELEROGRAM 


b.  RESPONSE  SPECTRUM,  6  PERCENT  DAMPING 

Figure  7.  Koyna,  India,  earthquake,  11  December  1967 
(transverse  component) 
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Figure  8.  Mode  shapes  for  the  San  Bernardino  tower 
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1 igure  9.  Mode  shapes  for  the  Briones  tower 
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Figure  10.  Response  of  the  San  Bernardino  tower  to  the  Taft, 

Calif.,  earthquake 
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igure  11.  Response  of  the  San  Bernardino  tower  to  the  Koyna , 

India,  earthquake 
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Figure  12.  Response  of  the  Briones  tower  to  the  Taft, 

Calif.,  earthquake 
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Figure  13. 

Response  of  the  Briones  tower 
India,  earthquake 

to  the  Koyna 
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where 

T  =  shearing  stress 

V^v  =  average  shear  on  a  cross  section 
A  =  cross-sectional  area 
CT  =  bending  stress 
M  =  moment 
r  =  average  radius 
I  =  moment  of  inertia 

As  these  thin-walled  approximations  are  only  valid  above  the  base  of 
the  structure,  the  results  are  shown  above  elevation  3254.75  for 
the  San  Bernardino  Tower  and  above  elevation  360.0  for  the  Briones 
Tower.  Such  sophisticated  models  may  be  required  for  the  final 
safety  evaluation  of  an  exceptionally  important  structure  subjected 
to  a  serious  seismic  hazard.  However,  from  the  cases  examined  with 
this  model,  two  important  conclusions  can  be  drawn  (Liaw  and  Chopra 
1973).  First,  the  effect  of  the  surrounding  water  is  approximately 
that  of  an  inertial  mass.  Second,  the  structure's  response  is  domi¬ 
nated  by  its  lower  flexural  inodes  of  vibration. 

Two-mode  added  mass  analysis 

10.  These  conclusions  suggest  that  the  seismic  response  of  an  in¬ 
take  tower  can  be  estimated  for  a  preliminary  analysis  by  considering 
the  first  two  flexural  modes  of  vibration  and  by  considering  an  appro¬ 
priate  amount  of  water  as  moving  with  the  structures.  Such  a  procedure 
has  been  developed  by  Chopra.-  In  this  procedure,  the  inclusion  of  two 
modes  has  been  found  to  be  necessary  to  adequately  predict  the  internal 
forces  induced  in  some  structures.  On  the  other  hand,  the  inclusion 


*  A.  K.  Chopra,  "Earthquake  Forces  for  Design  of  Intake-Outlet 

Towers,"  manuscript  submitted  to  the  U.  S.  Army  Engineer  Waterways 
Experiment  Station,  Vicksburg,  Miss.,  1981. 
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of  more  than  two  modes  has  been  judged  to  be  inappropriate  in  light  of 
the  approximation  of  hydrodynamic  interaction.  The  steps  in  the  pro¬ 
cedure  are  as  follows: 

±-  Define  the  structural  properties  of  the  tower-reservoir 
system: 

(1)  The  virtual  mass  per  unit  length  is: 


m(z)  =  m  (z)  +  m.(z)  +  m  (z) 

0X3 


where 

m  (z)  =  mass  of  the  tc.;er  by  itself 
o  J 

nu(z)  =  mass  of  the  water  inside  the  tower 

m  (z)  =  added  mass  to  represent  hydrodynamic  effects 
of  the  surrounding  water.  This  mass  can  be 
estimated  with  the  aid  of  Figure  14  which 
represents  the  effects  for  the  first  mode 
of  a  uniform  tower 

(2)  The  flexural  stiffness  per  unit  length  is  El(z)  . 

(3)  Damping  ratios  from  an  experimental  study  indicate 
that  5  percent  may  be  a  reasonable  value  for  this 
parameter  (Rea,  Liaw,  and  Chopra  1975). 

b.  Compute  the  periods,  T^  and  T^  ,  and  mode  shapes,  4» ^  Cz ) 
and  <J>2 Cz!'  ,  of  the  first  two  natural  modes  of  vibration. 
Rayleigh's  method  (Biggs  1964)  is  particularly  suited  for 
this  computation.  Alternately,  these  dynamic  character¬ 
istics  can  be  estimated  more  quickly  but  with  less  ac¬ 
curacy  from  Figures  15  and  16  in  which 

H 


is  the  average  virtual  mass  per  unit  length.  These 
design  aids  were  constructed  based  on  properties  of 
nonuniform  cantilever  beams  and  those  of  several 
actual  towers.  The  accuracy  of  these  aids  and  of 
Rayleigh's  metnod  for  the  San  Bernardino  and  Briones 
towers  can  be  seen  in  Figures  8  and  9,  respectively. 

c.  Compute  the  maximum  response  in  the  first  and  second 
modes  of  vibration: 


DISTANCE  ABOVE  BASE 


~  -1.0  \  0.75\  0.50 \  0.30  \0.20  \  0. 10\  0.05 


ADDED  MASS  RATIO 


w  2 

-  7r  ri 
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Figure  14.  Added  mass  representing  hydrodynamic  effects  of  surrounding 

water 

(1)  For  the  period  T  and  damping  ratio  4n  t  read 
the  pseudoabsolute  acceleration  S  from  the 
appropriate  design  responr  spectrum  as  discussed 
in  paragraph  6. 

(2)  Compute  the  equivalent  lateral  forces  from 


;  € 


I' 


L 

f  (z)  =  rp-  S  m(z)<j>  (z) 
n  M  an  Tn 

n 

where 

H 

L n  =  /  m(z)4>n(z)  dz 
o 

H 

M  -  /  m(z)4>^(z)  dz  ,  n  =  1,  2 
o  d 


c 


'a 


* 


(3)  Compute  the  internal  forces  (shears  and  moments)  at 
any  cross  section  by  a  static  analysis  of  the  tower 
subjected  to  the  equivalent  lateral  forces: 

H 

Vn(z)  =  f  fn(D  d(£) 
z 

H 

V2)  =  /  fn(£)(£  '  Z)  d(^  *  n  =  1,  2 

z 

\.here  z  <  £  <  H 

d.  Estimate  the  probable  maximum  shear  V(z)  and  moment 

M(z)  at  any  cross  section  by  combining  the  modal  maxima 
V  (z)  and  M  (z)  in  accordance  with  the  root  mean 
square  equations  of  Clough  and  Penzien  (1975): 

V(z)«  zT+V^T) 

M(z)~  -^Mj(z)  +  M*(z) 

11.  In  Figures  8-13,  the  results  of  analyzing  four  different 
tower-earthquake  permutations  by  this  t.wo-modal  approximation  procedure 
are  presented.  The  accuracy  of  the  procedure,  when  viewed  in  the  con¬ 
text  of  the  previously  discussed  refined  analysis,  may  be  adequate  for 
preliminary  design  evaluation.  However,  the  approximation  of  hydro- 
dynamic  effects  by  the  same  inertial  mass  in  both  modes  of  vibration 
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and  the  effects  of  tower  nonuniformity  should  be  examined  further 
before  this  observation  is  accepted  as  a  general  conclusion.  For 
towers  in  which  the  seismic  stresses  predicted  by  the  method  are  small 
compared  to  frequently  anticipated  stresses,  this  analysis  may  be  all 
that  is  required  for  final  design  evaluation  as  well.  The  steps  of 
this  approximate  procedure  for  a  particular  tower-earthquake  combina¬ 
tion  are  illustratively  detailed  in  Appendix  A. 

Montes-Rosenblueth  analysis 

12.  The  seismically  induced  shear  and  moment  distribution  for 
intake  towers  can  be  estimated  even  more  simply,  but  often  somewhat 
more  conservatively,  by  an  adaptation  (Chopra  1981)  of  a  procedure 
proposed  for  chimney  stacks  (Montes  and  Rosenblueth  1968).  The  simple 
expressions  for  these  distributions  were  developed  from  an  analysis  of 
a  uniform  flexural  cantilever  beam  excited  by  idealized  random  sup¬ 
port  motions.  The  procedure  starts  with  an  estimate  of  the  funda¬ 
mental  period  T^  of  the  intake  tower  including  the  mass  of  the  water 
inside  and  the  virtual  mass  of  the  surrounding  fluid.  Next,  shears 
and  moments  associated  with  a  flat  response  spectrum,  which  bounds 
the  design  response  spectrum  as  shown  in  Figure  17,  are  computed  from: 


where 

S  ~  max[S  (T) ] 

oi  3 

T  <  Tj 

H 

W  =  total  weight  =  g  J  m(z)  dz 

o 

Additionally,  shears  and  moments  associated  with  a  bounding  hyperbolic 
spectrum  having  a  low  period  cutoff  as  shown  in  Figure  17  are  calcu¬ 
lated  from 
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Since,  these  two  idealized  spectra  envelop  the  design  spectrum,  appro¬ 
priate  values  of  shear  and  moment  at  every  elevation  are  provided  by 
the  minimum  of  the  two  values  as  shown  for  the  values  of  shear  in 
Figure  18.  The  results  of  analyzing  four  different  tower-earthquake 
cases  by  this  procedure  are  shown  in  Figures  8-13.  The  procedure  can 
be  seen  to  be  quite  conservative  on  occasion  compared  to  the  previously 
discussed  two-mode  and  finite  element  methods.  However,  the  procedure 
is  considerably  simpler  to  apply  than  the  others  as  can  be  observed 


Figure  18.  Explanatory  sketch  to  obtain  the  Montes- 
Rosenblueth  (M-R)  estimate  of  shears 


from  the  sample  calculations  for  a  particular  tower-earthr jake  combina¬ 
tion  detailed  in  Appendix  B. 

Sequence  of  analyses 

13.  Thus,  in  an  evaluation  of  the  seismic  safety  of  a  tower  de¬ 
sign  or  an  existing  structure,  the  Montes-Rosenblueth  procedure  might 
be  used  for  a  first  analysis  since  ’t  is  simple  and  conservative .  If 
the  results  are  satisfactory,  th  ;  ana  lysis  need  not  be  carried  further. 
However,  if  the  results  are  unacceptable,  the  additional  work  required 
to  perform  the  more  accurate  two-mode  udded  mass  analysis  might  then  re 
warranted.  While  the  two-mode  added  mass  analysis  gives  shears  and 
moments  closer  to  those  obtained  by  the  refined  analysis  using  the 
EATSW  computer  program  than  the  Montes-Rosenblueth  method,  the  two-mode 
analysis  may  not  always  be  more  conservative  than  the  refined  analysis 
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using  EATS\.!  as  shown  in  Figures  10-13.  If  these  results  are  also  un¬ 
satisfactory,  a  refined  analysis  including  the  effects  of  hydrodynamic 
interaction,  such  as  that  of  Liaw  ard  Chopra  (1973),  might  be 
appropriate. 

Practical  Treatment  of  Asymmetry,  Foundation,  and  Embedment 

14.  The  analytical  methods  of  the  previous  section  assumed  the 
intake  tower  to  be  axisymmetric ,  with  a  rigidly  fixed  base,  and  free¬ 
standing.  In  practice,  considerable  departures  from  these  conditions 
occur.  Suggestions  for  the  treatment  of  these  departures  are  given  in 
the  following  paragraphs.  As  further  research  must  be  conducted  to 
rationally  defend  these  suggestions,  they  should  be  followed  with  con¬ 
siderable  caution. 

15.  For  reasons  of  operational  accessibility ,  some  intake  towers 
have  noncircular  cross  sections.  To  our  knowledge,  no  refined  analysis 
of  an  asymmetric  tower-reservoir  system  has  been  published.  However, 
in  one  case,*  a  rectangular  tower  was  analyzed  using  an  approximate 
axisymmetric  modal  method  with  a  slight  modification.  The  modifica¬ 
tion  replaced  each  rectangular  cross  section  with  an  equal  outer  area 
circular  one  to  calculate  only  the  added  masr;  representing  the  hydro- 
dynamic  interaction.  In  the  remainder  of  this  modal  analysis,  the 
properties  of  the  actual  rectangular  cross  section  were  used.  The 
appropriateness  of  this  modification  remains  to  be  evaluated  through 

a  more  refined  analysis. 

16.  Of  course,  no  intake  tower  base  is  absolutely  restrained, 
and  the  interaction  among  an  intake  tower,  its  reservoir,  and  its 
foundation  remains  a  viable  topic  for  future  research.  Practically 
speaking,  from  the  results  of  similar  research  for  dams  (Chopr3, 
Chakrabarti,  and  Gupta  1980),  it  is  expected  that,  if  the  foundation 


*  P.  F.  Mlakar,  "Preliminary  Seismic  Analysis  of  Cerillos  Tower," 

Letter  Report  submitted  to  the  U.  S.  Army  Engineer  District,  Jackson¬ 
ville,  from  the  U.  S.  Army  Engineer  Waterways  Experiment  Station, 
Vicksburg,  Miss  .  i'J/9. 
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stiffness  is  not  large  relative  to  the  tower  stiffness,  the  tower's 
natural  periods  will  be  longer  than  those  computed  by  the  fixed-base 
methods.  This  can  increase  or  decrease  the  dynamic  response  depend¬ 
ing  on  the  period  characteristics  of  the  site's  design  earthquake. 

This  effect  should  be  qualitatively  assessed  in  comparing  the  fixed- 
base  internal  forces  with  the  criteria  discussed  in  the  next  section. 

17.  Sometimes,  intake  towers  are  partially  .embedded  in  an  earth 
dam.  Recently,  the  effects  of  the  partial  embedment  of  a  specific 
structure  were  examined  with  a  rather  sophisticated  nonlinear  soil- 
structure  interaction  analytical  procedure.*  A  thorough  parametric 
investigation  of  the  seismic  behavior  of  such  structural  systems  is 
yet  to  be  performed.  However,  in  some  such  situations,  it  may  be 
possible  to  simply  and  conservatively  ignore  the  presence  of  the  sur¬ 
rounding  earth  and  ensure  the  seismic  safety  through  the  free-standing 
methods  of  the  previous  section.  If  these  results  are  unacceptable, 
the  free-standing  methods  might  be  modified  to  include  an  equivalent 
linearly  elastic  restraint  along  the  spatial  extent  of  partial 
embedment . 


Criteria 


18.  Once  the  distributions  of  seismically  induced  shear  and 
moment  have  been  determined  by  any  of  the  foregoing  preliminary  or  re¬ 
fined  procedures,  these  distributions  should  be  compared  to  the  capacity 
of  the  intake  tower  cress  sections.  In  this  comparison,  the  stresses 
from  the  static  loads  should  be  combined  with  the  stresses  from  the 
earthquake  loading.  The  loading  cases,  including  earthquake,  to  be  ex¬ 
amined  ire  contained  in  paragraph  3-07,  page  26,  of  Engineer  Manual 
(EM)  1110-2-2400  (Office,  Chief  of  Engineers  1964).  The  strength  de¬ 
sign  criteria  for  determining  the  structural  capacity  of  the  intake 
tower  are  contained  in  Engineer  Technical  Letter  (ETL)  1110-2-265 

*  Civil  Systems  Incorporated,  "Dynamic  Analysis  of  Structures  of 
Isabella  Dams,  Kern  County,  California,"  Report  to  the  U.  S.  Army 
Engineer  District,  Sacramento,  Calif. 
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(Office,  Chief  of  Engineers  1981).  An  example  of  the  strength  design 
of  an  intake  tower  is  contained  in  the  report  by  Liu  1980.  The  cri¬ 
teria  for  the  stability  of  the  tower  are  also  contained  in  para¬ 
graph  3-07  of  EM  1110-2-2400,  except  of  the  following:  the  seismic 
coefficients  for  the  overturning  and  sliding  stability  in  para¬ 
graph  3-07  have  been  superseded  by  those  listed  in  Engineer  Regula¬ 
tion  (ER)  1110-2-1806  (Office,  Chief  of  Engineers  1977)  and  the  cur¬ 
rent  sliding  analysis  criteria  are  contained  in  ETL  1110-2-256 
(Office,  Chief  of  Engineers  1983). 

19.  The  elastic  structural  response  implied  by  the  above  cri¬ 
teria  seems  appropriate  for  those  intake  towers  which  are  critically 
important;  i.e.,  in  those  cases  where  the  tower  would  be  needed  for 
a  controlled  release  of  the  reservoir  to  repair  any  seismic  damage  in 
the  damming  structure.  The  high  cost  associated  with  these  criteria 
for  a  noncritical  tower  loaded  by  a  credible  but  extremely  unlikely 
large  earthquake  suggests  the  possibility  of  inelastic  criteria  for 
such  structures.  However,  there  presently  exists  insufficient  in¬ 
formation  about  the  ductility  of  intake  towers  on  which  to  establish 
generally  applicable  inelastic  criteria.  Thus,  any  relaxation  of  the 
criteria  cited  in  the  preceding  paragraph  would  have  to  be  justified 
o«t  a  case-by-case  base. 


Additional  Research  Needed 


20.  As  mentioned  in  the  foregoing,  further  research  is  needed 
to  analyze  the  seismic  behavior  of  intake  towers.  Only  rather  modest 
efforts  would  be  required  to  examine: 

a.  The  effects  of  nonuniformity  on  the  added  mass  repre¬ 
sentation  of  hydrodynamic  interaction  in  Figure  14. 

b.  The  effect  of  representing  second-mode  hydrodynamic 
interaction  by  the  first-mode  added  mass  approximation 
in  Figure  14. 

c.  The  effect  of  assuming  contained  water  to  move  rigidly 
with  the  structure. 


d.  The  effect  of  approximating  modal  characteristics  with 
the  standard  shapes  in  Figures  15  and  16. 

e.  The  conditions  under  which  one,  two,  or  more  modes  of 
vibration  are  adequate  to  estimate  response. 

21.  More  extensive  studies  are  necessary  to  establish  the  ef¬ 
fects  of: 

a.  Asymmetry  on  hydrodynamic  interaction. 

b.  Three-dimensional  excitation  on  response. 

c.  Foundation  flexibility  on  behavior. 

d.  Structure  embedment  in  an  earth  embankment  on  response. 

e.  Inelastic  behavior  on  design  forces. 

Summary 


22.  Dynamic  methods  of  analysis  have  been  discussed  which  more 
rationally  represent  the  seismic  response  of  free-standing,  fixed-base 
intake  tower  reservoir  systems  than  the  seismic  coefficient  method. 

As  these  methods  explicitly  include  only  the  most  important  effects, 
they  should  be  used  cautiously  in  practice.  Further  investigations 
are  warranted  to  adequately  quantify  the  influence  of  simplifying 
assumptions  made  in  such  analyses. 
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Appendix  A:  Two-Mode  Added  Mass  Approximate  Analysis  of  the 
San  Bernardino  Tower  Excited  by  the  Taft,  Calif., 
Earthquake,  21  July  1952  (S69E  Component) 


1.  This  appendix  details  a  two-mode  added  mass  approximate 
analysis  of  the  San  Bernardino  tower  excited  by  the  S69E  component  of 
the  21  July  1952,  Taft,  Calif.,  earthquake. 

2.  The  material  properties  of  the  San  Bernardino  Tower  are 
assumed  as  follows: 

a.  Modulus  of  elasticity  of  concrete  =  4.5  million  psi. 

b.  Modulus  of  elasticity  of  steel  =  30  million  psi. 

c.  Poisson's  ratio  =  0.17. 

d.  Unit  weight  =  155  pcf. 

e.  Damping  =  5  percent  of  critical. 

3.  Figure  A1  and  Table  A1  show  the  discretization  of  this 

structure  for  the  computation  of  the  natural  periods,  Tj  ,  ,  and 

mode  shapes,  4^  ,  .  These  were  computed  using  a  Tektronix  4051 

BASIC  minicomputer  program  for  a  nonuniform  cantilever  beam  analysis. 
The  interactive  input-output  is  presented  on  pages  A4-A13,  while  the 
source  listing  is  given  on  pages  A14-A19. 


..  SAN  BERNARDINO  TOWER  b  RAYLEIGH  METHOD  GRID 


Figure  Al.  Discretization  of  the  San  Bernardino  tower  for  the 
Rayleigh's  method  calculation 


Table  A1 


Structural 


Properties  of  the  San  Bernardino  Tower 


Length 

in. 

Mass 

lb-sec^/in . 

Stiffness 
1015  lb-in 

102.00 

12,778.0 

55.50 

102.00 

8,752.0 

29.06 

102.00 

5,891 . 1 

13.52 

102.00 

3,537.7 

5.036 

396.72 

6,500.5 

1.073 

396.72 

6,410.7 

1.073 

409.56 

5,526.4 

1.073 

334  90 

1 ,980.8 

1.073 

149.16 

374.9 

0.50824 

149.04 

374.6 

0.50824 

2293.20 


52,126.7 


FIRST  AND  SECOND  FRENQUENCIES  OF  NON-UNIFORM 
CANTILEVER  USING  RAYLEIGH'S  METHOD 
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4.  Corresponding  to  the  periods 


T1  =  ?  =  =  °‘454  sec  T2  =  TO?  =  °'0888  sec 

the  spectral  accelerations 

S3l  =  0.432  g  Sa2  =  0.216  g 

were  then  obtained  from  Figure  6b  of  the  main  text. 

5.  Lateral  forces  f^(z)  acting  at  the  centroid  of  each  segment 
were  then  found  as  follows: 

Mx  =  Im(z)<t>^(z) 

2  2 
=  12,778(3.745  x  io"5)  +  8752(3.66  x  10-4) 

-3  2 

+  5891.1(1.262  x  io  °) 

2 

+  3537.7(3.339  x  10 ’3)  +  6500. 5(0. 0319)2  +  6410. 7(0. 185)2 

+  5526. 4(0. 426)2  +  1980 .8(0 . 7007)2  +  374.9(0.892)2 
+  374.6  (l)2 
Mj  =  2874.41 
Lj  =  ZaCzjf^z) 

=  5867.48 

Sa,  =  0.«2g  (386-4  f8/se-c2)  =  166.92  in. /sec2 

fj(z)  =  —  Sa1m(z)<t»1(z) 

M1 

S 

Sa.  =  340.73 
M1  1 

f  (51)  =  340. 73(1.2, 778) (3- 745  x  io"5  =  163.05  lb 
fj(153)  =  340 . 73(8752) (3.66  x  10'4)  =  1091.44  lb 
f  (255)  =  340. 73(5891. 1) (1 . 262  x  io"3)  =  2533.18  lb 


f j(357)  =  340.73(3537.7)(3.339  x  IO-3)  =  4024.84  lb 
^(606.36)  =  340. 73(6500. 5)(0. 0319)  =  70,655.95  lb 
^(1003.08)  =  340. 73(6410. 7)(0. 185)  =  404,099.64  lb 
^(1406.22)  =  340. 73(5526. 4)(0. 426)  =  802,164.05  lb 
f a (1803)  =  340. 73(1980.8)(0.7007)  =  472,916.02  lb 
fx(2069.58)  =  340. 73(374. 9)(0. 892)  =  113,944.03  lb 
^(2218.68)  =  340.73(374.6)0)  =  127,637.72  lb 
Lateral  forces  f2(z)  were  similarly  computed. 

M2  =  2m(z)<}>2(z) 

-4  2  o  2 

=  l2,778(-2.232  x  io  H)  +  8752(-2.l44  x  io""3) 

-3  2  9 

+  5891. 1 (- 7 .06  x  10  j  +  3 537 . 7 (-0 . 01737)2 
+  6500. 5 (-0.09347) 2  +  6410. 7(-0.3244)2  +  5526. 4(-0. 2409) 
+  1980. 8(0. 2576)2  +  374. 9(0 . 7189)2  +  374. 6(1)2 
M2  =  1753.34 

L2  =  2m(z)4>2(z)  =  -2988.83 
Sa2  =  0.216  B  (mA^idS«i)  =  83.46  i„./sec2 

h 

=  Sa2n(zH2(z) 

^2 

g-  Sa2  =  -142.27 

f 2  (5 1 )  =  (-142. 27)02, 778)  (-2. 232  x  io'4)  =  405.77  lb 
f 2 ( 153 )  =  (-142.27) (8752)(— 2. 144  x  IO-3)  =  2669.67  lb 
f2(255)  =  (-142.27) (5891 . 1) (-7.06  x  io"3)  =  5917.35  lb 
f 2 (357 )  =  (-142. 27)(3537. 7)(-0. 01737)  =  8742.72  lb 
0,(606.36)  =  (-142. 27) (6500.5) (-0 . 09347)  =  86,446.01  lb 


f2(1003.08)  =  (-142. 27) (6410 . 7) (-0 . 3244)  =  295,877.70  lb 
£2(14G6.22)  =  (-142. 27) (5526.4) (-0. 2409)  =  189,405.49  lb 
f 2 ( 1803 )  =  (-142. 27)(1980.8)(0. 2576)  =  -75,593.87  lb 
f2(2069.58)  =  (-142. 27) (374. 9) (0 . 7189)  =  -38,344.00  lb 
f2(2218.68)  =  (-142.27) (374. 61) (1)  =  -53,294.36  lb 

6.  Shears  V^(z)  were  then  calculated  from  internal  equilibrium. 


V^O)  =  163.05  +  1091.44  +  2533.18  +  4024.84  +  70,655.95  +  404,099.64 
+  802,164.05  +  472,916.02  +  113,944.03  +  127,637.72 
=  1,999,229.92  lb 
Vj (102)  =  1,999,066.87  lb 
Vj (204)  =  1,997,975.43  lb 
VjOOO)  =  1,995,442.25  lb 
Vj (408)  =  1,991,417.41  lb 
V1(804.72)  =  1,920,761.46  lb 
Vj (1201 . 44)  =  1,516,661.82  lb 
VjdOll)  =  714,497.77  lb 
VJU995)  =  241,581.75  lb 
Vx (,.144.16)  =  127,637.72  lb 
V.  (2293.2)  =  0 


Bending  moments  M^(z)  were  also  computed  from  equilibrium. 
(0)  =  163.05(51  -  0)  +  1091.44(153  -  0)  +  2533.18(255  -  0) 
+  4024.84(357  -  0)  +  70,655.95(606.36  -  0) 

+  404,099.64(1003.08  -  0)  +  802,164.05(1406.22  -  0) 

+  472,916.02(1803  -  0)  +  113,944.03(2069.58  -  0) 

+  127,637.72(2218.68  -  0)  =  2,950,135,600  in. -lb 


M1(102)  =  1091.44(153  -  102)  +  2533.18(255  -  102) 

+  4024.84(357  -  102)  +  70,655.95(606.36  -  102) 

+  404,099.64(1003.08  -  102)  +  802,164.05(1406.22  -  102) 
+  472,916.02(1803  -  102)  +  113,944.03(2069.58  -  102) 

+  127,637.72(2218.68  -  102)  =  2,785,574.064  in. -lb 
(204)  =  2533.18(255  -  204)  +  4024.84(357  -  204) 

+  70,655.95(606.36  -  204)  +  404,099.64(1003.08  -  204) 

+  802,164.05(1406.22  -  204)  +  472,916.02(1803  -  204) 

+  113,944.03(2069.58  -  204)  +  127,637.72(2218.68  ••  204) 
=  2,542,373,307  in. -lb 

1^(306)  =  4024.84(357  -  306)  +  70,655.95(606.36  -  306) 

+  404,099.64(1003.08  -  306)  +  802,164.05(1406.22  -  306) 
+  472,916.02(1803  -  306)  +  113,944.03(2069.58  -  306) 

+  127,637.72(2218  -  306)  =  2,338,708,971  in. -lb 
M1(408)  =  70,655.95(606.36  -  408)  +  404,099.64(1003.08  -  408) 

+  802,164.05(1406.22  -  408)  +  472,916.02(1803  -  408) 

+  113,944.03(2069.58  -  408)  +  127,637.72(2218.68  -  408) 
=  1,954,598,833  in. -lb 

Mj (804 . 72)  =  404,099.64(1003.08  ••  804.72)  +  802,164.05(1406.22 
-  804.72)  +  472,916.02(1803  -  804.72) 

+  113,944.03(2069.58  -  804.72) 

+  127,637.72(2218.68  -  804.72) 

=  1,359,359,362  in. -lb 
M1(1201.44)  =  802,164.05(1406.22  -  1201.44) 

+  472,916.02(1803  -  1201.44) 
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+  113,944.03(2069.58  -  1201.44)  +  127,637.72(2218.68 
-  1201.44) 

=  677,512,080  in. /lb 

M2 (1611 )  =  472,916.02(1803  -  1611)  +  113,944.03(2069.58  -  1611) 

+  127,637.72(2218.68  -  1611) 

=  220,615,219  in. -lb 

(1995)  -  113,944.03(2069.56  -  1995)  +  127,637.72(2218.68  -  1995) 

=  37,047,951  in-lb 

1^(2144.16)  =  127,637.72(2218.68  -  2144.16) 

=  9,511,563  in. -lb 
Mx(2293.2)  =  0  in. -lb 

The  shears  V2(z)  and  M2(z)  were  found  similarly. 

v2(0)  =  405.77  +  2669.67  +  5917.35  +  8742.72  +  86,446.01 
+  295,877 .70  +  189,405.49  -  75,593.87  -  38,344.00 
-  53,194.36 
=  422,232.48  lb 
V2(102)  =  421,826.71  lb 
V2(204)  =  419,157.04  lb 
V2(306)  =  413,239.69  lb 
V2(408)  =  404,496.97  lb 
V2(804.72)  =  318,050.96  lb 
V2(1201.44)  =  22,173,26  lb 
V2 (1611)  =  -167,232.23  lb 
V2(1995)  =  -91,638.36  lb 
V2(2144.16)  -  -53,294.4  lb 
V2(2293.2)  =  0  lb 


M2(102) 


M2(204) 


M2(0)  =  405.77(51  -  0)  +  2669.67(153  -  0)  +  5917.35(255  -  0) 

+  8742.72(357  -  0)  +  86,446.01(606.36  -  0) 

+  295,877.70(1003.08  -  0)  +  189,405.49(1406.22  -  0) 

-  75,593.87(1803  -  0)  -  38,344.00(2069.58  -  0) 

-  53,294.36(2218.68  -  0) 

=  286,716,569  in. -lb 

=  2669.67(153  -  102)  +  5917.35(255  -  102)  +  8742.72(357  -  102) 
+  86,446.01(606.36  -  107)  +  295,877.7(1003.08  -  102) 

+  189,405.49(1406.22  -  102)  -  75,593.87(1803  -  102) 

-  38,344.00(2069.58  -  192)  -  53,294.36(2218.68  -  102) 

=  243,669,551  in. -lb 

=  5917.35(255  -  204)  +  8742.72(357  -  204) 

+  86,446.01(606.36  -  204)  +  295,877.7(1003.08  -  204) 

+  189,405.49(1406.22  -  204)  -  75,593.87(1803  -  204) 

-  38,344.00(2069.58  -  204)  -  53,294.36(2218.68  -  204) 

=  200,779,380  in. -lb 

=  8742.72(357  -  306)  +  86,446.01(606.36  -  306) 

+  295,877.7(1003.08  -  306)  +  189,405.49(1406.22  -  306) 

-  75,593.87(1803  -  306)  -  38,344.00(2069.58  -  306) 

-  53,294.36(2218.68  -  306) 

=  158,327,146  in. -lb 

=  86,446.01(606.36  -  408)  +  295,877.7(1003.08  -  408) 

+  189,405.49(1406.22  -  408)  -  75,593.87(1803  -  408) 

-  38,344.00(2069.58  -  408)  -  53,294.36(2218.68  -  408) 

=  116,622,576  in. -lb 


M2(306) 


M2(408) 


O 
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(804.72)  -  295,877.7(1003.08  -  804.72)  <•  189,405.49  (1406.22  -  804.72) 

-  75,593.87(1803  -  804.72)  -  38,344.00(2069.58  -  804.72) 

-  5 j., 294. 36 (22 18. 68  -  804.72) 

-  -26,702,031  in. -lb 

M2 (1201 .44)  -  189,405.49(1406.22  -  1201.44)  -  75,593.87(1803 

-  1201.44)  -  38,344(2069.58  -  1201.44)  -  53,294.36(2218.68 

-  1201.44) 

-  -94,188,907  in. -lb 

M2 (1611)  -  -75,593.87(1803  -  1611)  -  38,344.00(2069.58  -  1611) 

-  53,294.36(2218.68  -  1611) 

-  -64,483,731  in. -lb 

M2 (1995)  -  -38,344.00(2069.58  -  1995)  -  53,294.36(2218.68  -  1995) 

-  -14,780,578  in. -lb 

M2(2144.16)  -  -53,294.36(2218.68  -  2144.16) 

=■  -3,971,496  in. -lb 
M2(2293.2)  -  0  in. -lb 

7.  Estimates  of  the  probable  maximum  shears  V(z)  were  then 
obtained  from  the  root  mean  square  equations. 

V(0)  *  -\J (1 ,999 ,229 .92) 2  +  (422, 732. 48)2  -  2,043,330.7  lb 
V(102)  -  2,043,087.4  lb 
V (204)  -  2,041,469.7  lb 
V (306)  -  2,037,782.3  lb 
V (408)  -  2,032,082.9  lb 
V (804 . 72)  -  1,946,915.8  lb 
V (1201.44)  -  1., 516, 823. 9  lb 


V ( 16 1 1 )  =  733,807.7  lb 
V(1995)  =  258,378,3  lb 
V (2144. 16)  =  138,317.3  lb 
V(2293. 2)  =  0 

Estimates  of  probable  maximum  moments  M(z)  were  similarly  computed. 

M(0)  = -y  2,950,135,600)2  +  (286 , 716 ,569)2  =  2,964,035,501  in. -lb 
M(  102)  =  2,796,211,315  in. -lb 
M(204)  =  2,550,289,080  in. -lb 
M(306)  =  2,344,062,102  in. -lb 
M(408)  =  1,958,074,928  in. -lb 
M(804, 72)  =  1,359,621,592  in. -lb 
M(1201 . 44)  =  684,027,900  in. -lb 
M(  161 1)  =  229,846,093  in. -lb 
M(1995)  =  39,887,544  in. -lb 
M(2144.16)  =  10,307,406  in. -lb 
M(2293 .2)  =  0 
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Appendix  B:  Montes-Rosenblueth  Approximate  Analysis  of  the 
San  Bernardino  Tower  Excited  by  the  Taft,  Calif., 
Earthquake,  21  July  1952  (S69E  Component! 


1.  Tj  =  0.454  seconds  (from  Appendix  A). 

2.  The  shears  and  moments  from  the  flat  spectrum  were  computed 


V(z)  =  0.647 


M(z)  =  0.461 


?)•['-  (i)! 


Using 


[sa(T)]  =  0.502  g 


T  <  Tj  =  0.454  sec 


(Figure  6b) 


H  =  2293.20  in. 

V(0)  =  6,541,921.8  lb 
V(  102)  =  6,541,346.1  lb 
V(204)  =  6,537,316.4  lb 
V(306)  '  6,526,378.5  lb 
V(408)  =  6,505,078.4  lb 
V(804 .72)  =  6,259,229.5  lb 
V ( 1 20 1 . 44)  =  5,601,176.7  lb 
V ( 1 6 1 1 )  =  4,273,801.0  lb 


20,141,756.9 

lb 

(Table  , 

(Table  , 

M(0)  = 

1.0689 

o 

o 

X 

in. 

-lb 

M( 102) 

=  9 

.984 

9 

x  10 

in. 

-lb 

M(204) 

=  9 

.295 

9 

x  10 

in. 

-lb 

M(306) 

=  8 

.623 

9 

x  10 

in. 

-lb 

M(408) 

=  7 

.967 

9 

x  10 

in . 

-lb 

M(804. 

72)  : 

=  5. 

5898  x 

109 

in. -lb 

M( 1201 . 

44)  = 

=  3. 

5113  x 

109 

in. -lb 

1 

'  1 


M( 161 1 )  =  1.7344  x  io9  in. -lb 


B1 


V C 1 995 )  =  2,234,590.8  lb 

V(2144. 16)  =  1,194,417.8  lb 

V(2293.2)  =  0 

3.  The  shears  and  moments 
computed  from 


M(z)  =  0.519 

V(0)  =  15,702,634  lb 
V(102)  =  15,509,818  lb 
V(204)  =  14,995,756  lb 
V(306)  =  14,249,975  lb 
V(408)  =  13,352,859  lb 
V(804 . 72)  =  9,612,385  lb 
V(1201 . 44)  =  7,269,238  lb 
V ( 16 1 1 )  =  6,888,661  lb 
V( 1995)  =  6,487,236  lb 
V(2144. 16)  =  5,206,164  lb 
V(2293.2)  =  0 


M C 1995 )  =  5.0124  x  10°  in. -lb 
M( 2 144 . 16)  =  1.771  x  108  in. -lb 
M(2293 .2)  =  0 

from  the  hyperbolic  spectrum  were 

-(b)2  - 6  25  (a) 2  ('  -  b)2. 

(¥) « (>  - 1) 

M(0)  =  1.2034  x  1010  in. -lb 
M(102)  =  1.1499  x  io10  in. -lb 
M(204)  =  1.0963  x  1010  in. -lb 
M(306)  -  1.0428  x  io10  in. -lb 
M(408)  =  9.893  x  109  in. -lb 
M(804 . 72)  =  7.811  x  io9  in. -lb 
M(1201 .44)  =  5.729  x  io9  in. -lb 
M(1611)  =  3.580  x  199  in. -lb 
M(1995)  =  1.5649  x  109  in. -lb 
M(2144 . 16)  =  7.8212  x  108  in. -lb 
M(2293 . 2)  =  0 


4.  The  shears  and  moments  from  the  flat  spectrum  should  be  used 
since  they  are  smaller  in  every  instance  than  those  from  the  hyperbolic 
spectrum. 


I 


In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a  facsimile  catalog 
card  in  Lihrary  of  Congress  MARC  format  is  reproduced 
below. 


Mlaker,  Paul  F. 

Seismic  analysis  of  intake  towers  /  by  Paul  F. 

Mlaker,  Patricia  S.  Jones  (Structures  Laboratory, 

U.S.  Army  Engineer  Waterways  Experiment  Station),  -- 
Vicksburg,  Miss.  :  The  Station  ;  Springfield,  Va,  ; 
available  from  NTIS,  1982. 

56  p.  in  various  pagings  ;  ill.  ;  27  cm.  -- 
(Technical  report  ;  SL.-82-8) 

Cover  title. 

"October  1982." 

Final  report. 

"Prepared  for  Office,  Chief  of  Engineers,  U.S.  Army." 
Bibliography:  p.  27. 

1.  Earthquake  prediction.  2.  Earthquakes  and  hvdraulic 
structures.  3.  Intakes  (Hydraulic  structures). 

4,  Seismologica'l  research.  5.  Seismology.  I.  Jones, 
Patricia  S.  II.  United  States.  Army.  Corps  of  Engineers. 


Mlaker,  Paul  F. 

Seismic  analysis  of  Intake  towers 


...  1982. 
(Card  2) 


Office  of  the  Chief  of  Engineers  III.  U.S.  Army 
Engineer  Waterways  Experin.  nt  Station.  Structures 
Laboratory.  IV.  Title  V.  Series:  Technical  report 
(U.S.  Army  Engineer  Waterways  Experiment  Station)  ; 
SL-82-8. 

TA7.W34  no .SL-82-8 


